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A new experiment is described to detect a permanent electric dipole moment of the proton with a 
sensitivity of 10~ 29 e-cm by using polarized “magic” momentum 0.7 GeV/c protons in an all-electric 
storage ring. Systematic errors relevant to the experiment are discussed and techniques to address 
them are presented. The measurement is sensitive to new physics beyond the Standard Model at 
the scale of 3000 TeV. 


INTRODUCTION 

One of the outstanding problems in contemporary el¬ 
ementary particle physics and cosmology is finding an 
explanation for the observed matter-antimatter asym¬ 
metry of our universe, known as baryogenesis. Within 
the framework of the Big Bang, it appears that a much 
greater degree of CP-violation than provided by the Stan¬ 
dard Model of particle physics is required. That suggests 
the necessary existence of New Physics (NP), with large 
CP-violating interactions as a key ingredient in under¬ 
standing the early universe. Identifying that NP source, 
central to our very existence, would be a major intellec¬ 
tual achievement. 

In the search for laboratory manifestations of new CP- 
violating effects, electric dipole moments (EDMs) play 
a crucial role. A non-zero permanent particle electric 
dipole moment (EDM) separately violates parity (P) and 
time reversal symmetry (T) [Tj. So, assuming CPT in¬ 
variance, CP must also be violated. Since Standard 
Model EDM predictions are much smaller than current 
experimental sensitivities, an observation of any parti¬ 
cle’s EDM with today’s technology would signal discov¬ 
ery of NP. If of sufficient strength, such a source could 
provide a possible explanation for baryogenesis. Some 
theories, which suggest that an EDM may be within 
experimental reach, include supersymmetry (SUSY) [5], 
left-right symmetry [3j, and multi-Higgs scenarios [3]. 
Here, we explore the possibility of a storage ring search 
and study of the proton EDM (d p ) at the unprece¬ 
dented level of 10 _29 e • cm, an advance by nearly 5 or¬ 
ders of magnitude beyond the current indirect bound of 
\d p \ < 7.9 x 10 _ 25 e-cm obtained using Hg atoms [5|. Ob¬ 


serving the EDM from different simple systems is neces¬ 
sary to identify the source of any NP [Bj . 

This dedicated direct proton EDM study at the level 
of 10 _29 e • cm is sensitive to a generic NP mass scale 
Anp with CP-violating phase <^>np roughly satisfying [7] 
(3000 TeV/ANp ) 2 tan (^np) > 1- For a phase of the 
order of 45 degrees, the 3000 TeV NP scale is being 
probed, while for NP generically parametrized by a scale 
of order 1 TeV, a relative phase sensitivity as small as 
</> NP ss 10“' would be reached. Many specific exam¬ 
ples of the outstanding probing power of a proton EDM 
study at the 10 _29 e • cm level exist. Here, we point out 
that it constrains the 0qcd parameter at 10 -13 , 3 orders 
of magnitude below the current neutron EDM bounds. 
A more timely illustration in the current LHC era is a 
potential CP-violating chiral phase induced by a loop 
induced Higgs to 2-photon coupling (the relative pseu¬ 
doscalar to scalar amplitudes, a measure of potential 
Higgs CP-violation). Such an effective coupling would 
lead to fermion EDMs via quantum loops, making for an 
overall effect of 2-loop order. The proton EDM experi¬ 
mental program envisioned in this paper would be sen¬ 
sitive to a relative pseudoscalar coupling of order 10 -3 , 
about 2 orders of magnitude below current electron EDM 
constraints. It should be noted that EDMs of the elec¬ 
tron, proton and neutron may actually represent our only 
practical access to the very small Higgs coupling to light 
(first generation) fermions, requiring of course the caveat 
that CP is violated at an observable level in the H 77 
interaction. 

Searching for a non-zero proton EDM in a dedicated 
storage ring presents an experimental opportunity to im¬ 
prove the current sensitivity by more than three orders 
of magnitude compared to the current neutron EDM ex- 
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perimental limit 0 . The method we describe is based on 
the frozen spin method and uses an all-electric lattice, di¬ 
rectly measuring spin precession due to a non-zero EDM 
in an electric field. The Storage Ring EDM Collabora¬ 
tion has made significant progress in developing the ex¬ 
perimental design for an all-electric EDM measurement. 
In this paper, we describe the fundamental experimental 
techniques and the specifications of the all-electric stor¬ 
age ring. We also present the systematic errors and the 
methods developed to address them. 


EXPERIMENTAL METHOD 


The EDM and magnetic moment in terms of the rest 
frame spin s are d = (7ye/2mc)s, and fi = (ge/2m)s 1 
respectively, where these relations define g and 77 ; G = 
(g — 2)/2 defined; and m is the particle mass. At rest, 
the spin precession of a particle in electric and magnetic 
fields E and B is governed by: 

ds „ , _ 

— = fi x B + d x E. (1) 

at 

For a particle with velocity (3 = v/c, relativistic factor 
7 = (l — v 2 /c 2 ) 1 ^ 2 , (3 ■ E = 0 and (3 ■ B = 0, the spin 
precesses relative to the momentum with angular velocity 
u) a + uj e PflT| , where: 


TO 


GB — [ G — 


7 2 — 1 


f3 x E 




( 2 ) 


Spin precession in a storage ring has been successfully 
used to establish a limit on a muon EDM nans]. For cer¬ 
tain E and B fields, the (g — 2) precession of the particle 
uj a vanishes mna. Thus, aside from an EDM contri¬ 
bution, the spin is frozen along the momentum direction. 
With B = 0 and the “magic” 7 = -y/l + 1 /G chosen, 
u3 a = 0 and cu e = (ge/2mc)B. The parameters for this 
condition in the case of a proton are shown in Table |T] 
In an all-electric storage ring chi eg, a radial electric 
field causes the spin to precess out of the storage plane 
linearly on the time-scale of the fill. Details of a storage 
ring proton EDM experiment are given in [20] . For long- 
lived, polarized beams, the gain in sensitivity by using 
the frozen spin method over indirect methods in CHI is 
several orders of magnitude [50] . 


TABLE I. “Magic” proton parameters to cancel the (g — 2) 
precession in an all-electric ring. 


g„ EH 

7 

0 

P 

E 

1.792847 

1.248107 

0.598379 

0.7007 GeV/c 

1.171 GeV 


Ring and beam parameters. The all-electric ring ge¬ 
ometry will include 40 sections of concentric cylindrical 
deflectors of 52.3 m bending radius, with 36 straight sec¬ 
tions of 2.7 m length and four straight sections of 20.8 m 
length, adding up to a 500 m circumference. The 2.7 m 
straight sections will include superconducting quantum 
interference devices (SQUIDS) as magnetometers and 
electrostatic alternating gradient quadrupoles, with two 
polarimeters placed in the longer straight sections. In¬ 
jection of the beams in opposite directions around the 
ring will occur in the remaining two 20.8 m straight sec¬ 
tions. The deflector electric field will be about 8 MV/m 
radially inward in the 3 cm spacing between the deflector 
plates. In principle it is possible to modify the shape of 
the deflector plates to include vertical focusing. The ver¬ 
tical bending specifications are strict, on the micrometer 
scale, and are under study. An electrostatic storage ring 
of this size would be more than ten times larger than any 
previous electrostatic ring [221124] . A simplified deflector 
lattice is shown in Figure |T| 

Experimental techniques. One hundred bunches of 
2.5 x 10 s vertically polarized protons will be injected in 
the clockwise (CW) direction and a similar number in the 
counter-clockwise (CCW) direction, circulating simulta¬ 
neously. The beams will be allowed to de-bunch and then 
re-bunch at the required frequency using a high-harmonic 
(h = 100 in this case) RF system. Using an RF solenoid, 
the spins of the protons will be rotated in the longitu¬ 
dinal direction, producing protons of both positive and 
negative helicities. The vertical polarization difference 
between early and late times will determine the average 
vertical spin precession rate. The considered ring and 
beam parameters are summarized in Table ^ 

The horizontal spin coherence time (SCT) of the stored 
beam is the time required for the RMS spread in spin an¬ 
gles to become one radian. Causes of a finite SCT include 
horizontal and vertical oscillations as well as longitudinal 
(energy) oscillations. An RF-cavity is required to keep 
the SCT much longer than a few ms due to momentum 
spread of the beam. The vertical component of the pro¬ 
ton spin grows linearly with time, at a rate on the order 
of nanoradians per second for an EDM detectable by the 
experiment. In practice, the linear growth of the vertical 
spin is limited by the SCT of the stored beam. Protons 
will be stored on the order of 10 3 s, yielding an early-to- 
late change of the vertical spin component on the order 
of microradians. A vacuum of better than 10 -10 Torr is 
required to keep the beam stored in the ring for 10 3 s. 

The electric field of 8 MV/m between the cylindri¬ 
cal deflectors is comparable to previous work which has 
achieved similar field strengths |25H27j . The effect of the 
fringe electric fields of the cylindrical deflectors has been 
investigated both analytically and by precision particle 
tracking [2S]. The bending radius of the plates has been 
adjusted to account for additional deflection due to fringe 
electric fields in the straight sections. 
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Intrabeam scattering (IBS) effects, increasing the 
stored beam phase-space parameters, set the time scale 
of the fill. The necessity of having sufficiently small IBS 
results in a moderate beam current limit. The resulting 
space-charge tune shifts are small. The beam-beam scat¬ 
tering effects are smaller and do not present a problem. 



FIG. 1. Top: An overview of an idealized ring with 16 cylin¬ 
drical deflector segments and straight sections. The deflector 
bending radius is Ro = 52.3 m and plate spacing is 3 cm. The 
electric field is directed radially inward between the plates. 
The spin and momentum vectors are kept aligned for the du¬ 
ration of storage. Bottom: A realistic lattice will include 40 
bending sections separated by 36 straight sections 2.7 m long 
each, with electrostatic quadrupoles in an alternating gradi¬ 
ent configuration, and four 20.8 m long straight sections for 
polarimetry and beam injection. It will also include SQUID- 
based magnetometers distributed around the ring, whose total 
circumference is 500 m. 


TABLE II. Ring and beam parameters of the proton EDM 
experiment. Proton beam parameters refer to each storage 
direction. 


Bending radius, Rq 

52.3 m 

Electrode spacing, d 

3 cm 

Electrode height 

20 cm 

Deflector shape 

cylindrical 

Radial E-field, Eq 

8 MV/m 

Number of straight sections 

40 

Straight section lengths 

2.7389 m, 20.834 m 

Polarimeter sections 

2 

Injection sections 

2 

SQUID-based magnetometer sections 

36 

Total circumference, C 

500 m 

Harmonic number h, RF frequency 

100, 35.878 MHz 

RF voltage, synchrotron tune Q s 

6 kV, 0.0066 

Particles per bunch 

2.5 Xl0 8 

Maximum momentum spread, (dp/p ) ma x 

4.6 X 10 ~ 4 

Horizontal beta function, /3 X)max 

47 m 

Vertical beta function, /3y, max 

216 m 

Horizontal dispersion function, D X)max 

29.5 m 

Horizontal tune, Q x 

2.42 

Vertical tune, Q y 

0.44 

Vertical emittance, evmax 

17 mm mrad 

Horizontal emittance, £Hmax 

3.2 mm mrad 

Slip-factor, 77 = a — I/ 7 2 

-0.192 


Polarimeters will be used to measure small changes in 
the vertical component of the beam polarization by scat¬ 
tering particles from a 6 cm thick carbon target. The 
stored beam will be led to collide with the target us¬ 
ing a number of possible alternative methods, e.g., by 
slowly lowering the vertical focusing strength, using a 
resonant slow extraction vertically, etc. About 99% of 
the time the beam particles undergo Coulomb scatter¬ 
ing, lose enough energy and leave the ring. Roughly 1% 
of the time the protons undergo spin-dependent nuclear 
elastic scattering, ending up on a detector located about 
a meter beyond the target. The scattering of the parti¬ 
cles in the up and down (left and right) directions will 
provide information on the horizontal (vertical) plane po¬ 
larization component. Detectors must be able to respond 
to charged-particle events with minimal dead time and 
small systematic errors. Types under consideration in¬ 
clude multi-resistive plate chambers, micro-megas cham¬ 
bers, gas electron multiplier chambers, and silicon detec¬ 
tors. A polarimeter design under consideration is shown 
in Figure [2j 
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FIG. 2. A possible layout of one-half of a proton EDM ring 
polarimeter, with the scale in centimeters. The carbon extrac¬ 
tion target is located on a motor drive. The beam is denoted 
by a short-dashed line. 

A precision tracking program using numerical inte¬ 
gration methods such as fourth-order Runge-Kutta and 
the Predictor-Corrector method m is used to inte¬ 
grate the beam and spin dynamics differential equa¬ 
tions. The Runge-Kutta and Predictor-Corrector inte¬ 
gration methods used in tracking are slow but accurate; 
they are kept at a step size on the order of 1-10 ps. 
The programs have been benchmarked to high accuracy 
against analytical estimates. Additional programs, such 
as UAL/ETEAPOT |311 152 ] or using the Hamiltonian 
approach [331 . which account for the momentum change 
of the particle due to motion in E-fields, have been devel¬ 
oped. Results from tracking simulations [SU [35] confirm 
that a particle distribution with (dp/p) = 2 x 10~ 4 will 
have a SCT of the order 10 3 s in an idealized lattice. 
The SCT of the stored beam was found to increase with 
increasing ring radius. 

For a uniform beam extraction rate, the statistical er¬ 
ror <Jd of the measurement has been both calculated an¬ 
alytically and confirmed numerically with Monte Carlo 
simulations to be: 


P IVtot, cltot/TSCT 

Taking parameter values |39] of the beam polarization 
P = 0.8, the analyzer power A = 0.6, E 0 = 8 MV/m over 
65% of the ring, N totiC = 5 x 10 10 particles per storage 
cycle, a particle detection efficiency of / = 0.0055, a 
total running time Ttot = 10 7 s per year, and a SCT of 
TgcT = 10 3 s, the statistical error for one year is 4 x 
l(U 29 e-cm. The polarimeter analyzing power peaks near 
the proton “magic” momentum, a fortuitous coincidence. 
The statistical error grows in proportion to the storage 
ring radius. That error will be reduced by modulating 
the data-taking rate, taking most of the data at early 
and late times. By further optimizing the ring lattice, 
reducing IBS and increasing SCT, we expect to be able 
to achieve a statistical sensitivity of < 2 x 10 _29 e-cm per 
year. 


SYSTEMATIC ERRORS 


Analytical estimates in combination with precision 
tracking allow the size of potential systematic errors to 
be estimated in several ways. Magnetic shielding, beam 
position monitors (SQUID-based BPMs sensitive to B- 
fields, plus button-BPMs sensitive to E-fields), and lat¬ 
tice alignment to better than 0.1 mm around the ring 
are sufficient to address the main systematic errors, i.e., 
radial B-fields and geometric phases [36]. 

Radial B-fields. The ring will be shielded from the 
Earth’s magnetic field as well as from noise through pas¬ 
sive shielding and feedback mechanisms. The presence of 
a net radial B-field, ( B r ), could mimic an EDM signal, 
producing a vertical spin precession 


eg (B r ) 

2 mj 2 


(4) 


An average radial field of 10 aT will cause a spin preces¬ 
sion at the sensitivity of the EDM experiment. However, 
the radial magnetic field would split the CW and CCW 
beams vertically. The vertical beam position Sy as a 
function of the modes of the radial B-field is: 


& y= Y. 

N =0 


(dcRoB r N 

E 0 (Q 2 y ~ N*) 


cos(N9 + pn), 


(5) 


where Q y is the vertical tune, which will oscillate between 
0.44 and 0.48 at a frequency on the order of 10 kHz. 

As noted above, the counter-rotating beams will be 
split vertically, with maximum separation equal to twice 
that given by Equation [5] For a 10 aT magnetic field, 
the split between the beams is on the order of a pi- 
cometer. Beam position monitors (BPMs) are required 
that can determine the vertical positions of the CW and 
CCW beams with picometer-scale resolution over the 
10' s duration of the experiment. SQUID magnetome¬ 
ters are suitable to measure the magnetic field resulting 
from the splitting of the CW and CCW beams. Low- 
temperature DC SQUIDs have demonstrated sensitivities 
down to 1 fT/-\/Hz [37, ,3BJ ■ Measuring an average split¬ 
ting between the CW and CCW beams on the picometer 
level is feasible with appropriate SQUID placement in the 
straight sections of the lattice. In principle, no magnetic 
shielding will be required with continuous BPM measure¬ 
ment around the ring. The finite number of detectors 
limits the magnetic field modes to which the system is 
sensitive by the Nyquist sampling theorem. The critical 
parameter is the detected average radial B-field and that 
can be wrong only when the mode is equal to the num¬ 
ber of the BPM locations around the ring as well as its 
integer multiples. However, it has been shown analyti¬ 
cally as well as by beam/spin tracking simulations that 
the SQUID-based BPMs are quite insensitive to higher 
harmonics of the radial B-fields. SQUIDS are sensitive 
to time dependent B-fields only. The ratio of the time 
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oscillating component of N = 0 to higher harmonics can 
be shown to be ( Q y /N ) 4 , making the high N harmonic 
contributions negligible. The shielding requirements are 
dominated by the so-called geometric phase effect (see 
below). For this reason, the total magnetic field will be 
shielded to below 10-100 nT at all points in the ring. 
Such shielding of the B-field is within the present state 
of the art [40]. 

SQUIDs in the straight-sections of the ring should be 
sufficient to measure the vertical splitting. A schematic 
of one such possible SQUID BPM station is shown in 
Figure [3] The vertical spin precession rate as a function 
of the detected radial B-field will be plotted, with the 
EDM signal corresponding to the DC offset in the vertical 
axis. 


Aluminum 
Mu-metal 
Dewar 
4.2 K 
Nb tube 
Pickup coil 
Beam 
Beam pipe 
RF screen 


FIG. 3. A schematic of a possible SQUID BPM station. The 
system is shielded with a superconducting Nb tube, A1 tube 
for RF-shield, and several mu-metal layers. 



Geometric phases. In three dimensions, spin rota¬ 
tions about different axes do not commute. This fact 
contributes to geometric phase-induced false EDMs, a 
significant systematic error in neutron EDM experi¬ 
ments El [Ml HI]. In a storage ring, geometric phases 
of spin dynamics may be cancelled more simply than 
in a neutron trap. Consider, for example, the iV-modes 
of two non-commuting perturbations: spin rotation fre¬ 
quency around the vertical axis, (Suv)n cos (Noj c t), 
and spin rotation frequency around the longitudinal 
axis, (5ujl)n cos (Noj c t + </>). The presence of these 
perturbations in the lattice leads [|5] to a spin rota¬ 
tion around the radial axis with a frequency Qr = 
\(Su;v)n{Scjl)n sincj)/(2Nuj c )\, thus imitating an EDM 
signal. It follows from this formula that in order to can¬ 
cel the geometric phase we only need to find, exper¬ 
imentally, a counter-perturbation for either of the two 


iV-modes of the perturbed fields. The same approach 
can be used for other modes. 

Perturbed E and B-fields induce distortions of the 
closed orbits, the detection of which will be used to check 
that the errors are being kept at an acceptable level [50] ■ 
With respect to achieving this level, analytic investiga¬ 
tion has shown that the geometric effect of B-fields split¬ 
ting the counter-rotating beams can be kept below the 
experimental sensitivity if the maximum B-field is kept 
below the 10-100 nT level everywhere; that E-field er¬ 
rors due to plate misalignments displacing both counter¬ 
rotating beams from their ideal location can be addressed 
by placing BPMs within 0.1 mm of the ideal orbit; and 
that small changes to deflector geometry will adequately 
compensate for systematic errors due to deflector fringe 

fields [58] ■ 

Polarimetry. An experiment investigating the manage¬ 
ment of geometric and rate-induced systematic errors in 
polarimetry was conducted at COSY-Jiilich [55] with a 
1.7 cm carbon block target. Large systematic errors con¬ 
sisting of position and angle changes to the beam were 
made deliberately in order to generate easily measurable 
effects. At the level of geometric and rate errors ex¬ 
pected for the proton EDM experiment, the results of 
the COSY-Jiilich study indicate how to reduce system¬ 
atic uncertainties in polarimetry to well below the level 
of sensitivity. Positive and negative helicity bunches will 
be stored in the same direction, and a combination of 
observables will be used to identify systematic errors due 
to non-linearities [M] ■ 

Vertical forces. The interaction of the CW and CCW 
beams may lead to a systematic error: if the two counter¬ 
rotating beams do not overlap completely, on average 
they will feel a vertical force from one another. The 
problem can be addressed as long as the SQUID BPMs 
are sensitive to the beam separation size, and feedback 
can be used to eliminate the signal. Any forces on the 
beams due to image charges on the top and bottom of 
the vacuum chamber will be minimized by using verti¬ 
cal metallic plates for almost the entire azimuthal extent 
of the ring. Results from numerical simulations indicate 
that the aspect ratio of the quadrupole plates can be cho¬ 
sen to reduce the effect when the counter-rotating beam 
intensities do not cancel exactly. Sextupole electric fields 
combined with different CW and CCW beam sizes can 
create a vertical splitting of the counter-rotating beams, 
setting the specifications for both. Introducing short runs 
in between regular-length runs will address systematic er¬ 
rors related to the vertical spin component of the beams 
being correlated to the protons’ phase space parameters. 

Other effects. There are potential systematic errors 
due to the gravitational field and rotation of the Earth. 
For example, there is a false-EDM signal due to the verti¬ 
cal E-fields being balanced by the force of gravity at our 
level of sensitivity. Taking the difference between sig¬ 
nals of the CW and CCW beams will cancel this effect. 
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Also, Coriolis and Sagnac [35] effects due to the rotation 
of the Earth have been found to be below the exper¬ 
imental sensitivity. The RF cavity will account for the 
slightly different travel times of the CW and CCW beams 
around the ring by equalizing the frequencies. However, 
the Sagnac effect may place an upper limit (more than 
10 3 s) on how long counter-rotating beams will be stored 
with longitudinal polarization. 

Spin resonances may also contribute to a false EDM 
signal, although this contribution is decreased in the 
frozen spin ring by a factor 1 /tsct- In addition, spin 
and beam resonances coincide in the frozen spin ring, so 
they will be dealt with together (at a later stage of the 
project). 

A summary of the main systematic errors in the ex¬ 
periment is given in Table |III| 


TABLE III. Main systematic errors of the experiment and 


their remediation. 


Effect 

Remediation 

Radial B-field 

SQUID BPMs with 1 fT/\/Hz sensitivity 
eliminate it. 

Geometric phase 

Plate alignment to better than 100 /rm, 
plus CW and CCW storage. Reducing 
B-field everywhere to below 10-100 nT. 
BPM to 100 fim to control the effect. 

Non-Radial E-field 

CW and CCW beams cancel the effect. 

Vert. Quad misalignment 

BPM measurement sensitive to vertical 
beam oscillation common to CW and 
CCW beams. 

Polarimetry 

Using positive and negative helicity pro¬ 
tons in both the CW and CCW directions 
cancels the errors. 

Image charges 

Using vertical metallic plates except in 
the quad region. Quad plates’ aspect ra¬ 
tio reduces the effect. 

RF cavity misalignment 

Limiting longitudinal impedance to 10kQ 
to control the effect of a vertical angu¬ 
lar misalignment. CW and CCW beams 
cancel the effect of a vertically misplaced 
cavity. 


CONCLUSIONS 

The Storage Ring EDM Collaboration has designed 
an experiment using the frozen spin method and a ded¬ 
icated storage ring to measure the proton EDM with 
an unprecedented sensitivity of lCP 29 e-cm. We are cur¬ 
rently developing prototypes to optimize the critical sys¬ 
tems of the experiment, which include magnetic shield¬ 
ing, SQUID-based BPMs, polarimeter and electric field 
plates. In parallel we are developing software for high 
precision and high efficiency spin and beam dynamics 
tracking. The proton EDM measurement will provide 
a valuable probe of new physics beyond the Standard 
Model. 
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